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Abstract 

We give new formulas on the total number of born particles in the stable birth- 
and-assassination process, and prove that it has a heavy-tailed distribution. We 
also establish that this process is a scaling limit of a process of rumor scotching 
in a network, and is related to a predator-prey dynamics. 
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1 Introduction 

Birth-and-assassination process 

The birth-and-assassination process was introduced by Aldous and Krebs [2], it is a 
variant of the branching process. The original motivation of the authors was then to 
analyze a scaling limit of a queueing process with blocking which appeared in database 
processing, see Tsitsiklis, Papadimitriou and Humblet [13] . In this paper, we show 
that the birth-and-assassination process exhibits some heavy-tailed distribution. For 
general references on heavy-tail distribution in queueing processes, see for example 
Mitzenmacher [9] or Resnick [12]. In this paper, we will not discuss this application. 
Instead, we will show that the birth-and-assassination process is also the scaling limit 
of a rumor spreading model which is motivated by network epidemics and dynamic 
data dissemination (see for example, [TO], [3J, [TT]). 

We now reproduce the formal definition of the birth-and-assassination process from 
[2J. Let N J = U£L N fc be the set of finite k-tuples of positive integers (with iV° = 0). Let 
{$„}, n G W, be a family of independent Poisson processes with common arrival rate 
A. Let {K n },n e W, be a family of independent, identically distributed (iid), strictly 
positive random variables. Suppose the families {$ n } and {K n } are independent. 
The particle system starts at time with only the ancestor particle, indexed by 0. 
This particle produces offspring at the arrival times of $0, which enter the system 
with indices (1), (2), ■ • • according to their birth order. Each new particle n entering 
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Figure 1: Illustration of the birth- and- assassination process, living particles are in red, 
dead particles in blue, particles at risk are encircled. 



the system immediately begins producing offspring at the arrival times of {$ n }, the 
offspring of n are indexed (n, 1), (n, 2), • • • also according to birth order. The ancestor 
particle is at risk at time 0. It continues to produce offspring until time D% = K$, 
when it dies. Let k > and let n = (n\, • • • ,nk-x,n^), n' = (n 1; nk-i). When a 
particle n' dies (at time D n i), n then becomes at risk; it continues to produce offspring 
until time D n = D n * + K n , when it dies. We will say that the birth-and-assassination 
process is stable if with probability 1 there exists some time t < oo with no living 
particle. The process is unstable if it is not stable. Aldous and Krebs [2] proved the 
following: 

Theorem 1 (Aldous and Krebs) Consider a birth-and-assassination process with 
offspring rate A whose killing distribution has moment generating function <ft. Suppose 
cf) is finite in some neighborhood ofO. If min u>0 Au _1 0(m) < 1 then the process is stable. 
If min u> o Am _1 0(m) > I then the process is unstable. 

The birth-and-assassination process is a variant the classical branching process. 
Indeed, if instead the particle n is at risk not when its parent dies but when the particle 
n was born, then we obtain a well-studied type of branching process, refer to Athreya 
and Ney [5]. The populations in successive generations behave as the simple Galton- 
Walton branching process with mean offspring equal to AEA0, and so the process is 
stable if this mean is less than 1. The birth-and-assassination process is a variation in 
which the 'clock' which counts down the time until a particle's death does not start 
ticking until the particle's parent dies. 

In this paper, we will pay attention to the special case where the killing distribution 
is an exponential distribution with intensity fi. By a straightforward scaling argument, 
a birth-and-assassination process with intensities (A, /x) and a birth-and-assassination 
process with intensities (A// -1 , 1) where the time is accelerated by a factor /1 have the 
same distribution. Therefore, without loss of generality, from now on, we will consider 
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B, a birth-and-assassination process with intensities (A, 1). As a corollary of Theorem 
[U we get 

Corollary 1 (Aldous and Krebs) If < A < 1/4, the process B is stable. If A > 
1/4, the process B is unstable. 

In the first part of this paper, we study the behavior of the process B in the stable 
regime, especially as A get close to 1/4. We introduce a family of probability measures 
{Pa}, A > 0, on our underlying probability space such that under P x , B is a birth-and- 
assassination process with intensities (A, 1). Let A G (0, 1/4), we define N as the total 
number of born particles in B (including the ancestor particle) and 

7(A) = sup {u > : E X N U < oo} . 

In particular, if < 7(A) < 00, from Markov Inequality, for all < e < 7(A), there 
exists a constant C > 1 such that for all t > 1, 

P\{N >t)< cr 7(A)+e . 

The number 7 may thus be interpreted as a power tail exponent. There is a simple 
expression for 7. 

Theorem 2 For all A G (0, 1/4), 

... l + VT^4\ 
7(A) 



1 - VI -4A 

This result contrasts with the behavior of the classical branching process, where for 
all A < 1: there exists a constant c > such that EAexp(ciV) < 00. This heavy tail 
behavior of the birth-and-assassination process is thus a striking feature of this process. 
Near criticality, as A j 1/4, we get 7(A) ~ 1, whereas as A J. 0, we find 7(A) ~ (2A) _1 . 
By recursion, we will also compute the moments of N. 

Theorem 3 (i) For all p > 2, E X N P < 00 if and only if A G (0,p(p+ 1)~ 2 ). 



(ii) //AG (0,1/4], 



E X N = (1) 

l + v / T^4A V ; 



E X N 2 = - . (2) 



(111) //AG (0,2/9), 

Theorem [3]^ i) is consistent with Theorem [2j A G {0,p(p + 1)~ 2 ) is equivalent to p G 
[1, (1 + Vl — 4A)(1 — Vl — 4A) _1 ) . Theorem [3](ii) implies a surprising discontinuity 
of the function A 1— > E X N at the critical intensity A = 1/4: lmTvri/4 E X N = 2. Again, 
this discontinuity contrasts with what happens in a standard Galton- Watson process 
near criticality, where for < A < 1, E X N = (1 — A) -1 . We will prove also that this 
discontinuity is specific to A = 1/4 and for all p > 2, lim A | p ( p+1 )-2 E A [iV p ] = 00. We 
will explain a method to compute all integers moments of N by recursion. The third 
moment has already a complicated expression (see §2.5.11) . From Theorem Etli), we 
may fill the gap in Corollary [TJ 
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Corollary 2 If A = 1/4, the process B is stable. 

In Section [2j we will prove Theorems [2] and [3] by exhibiting a Recursive Distribu- 
tional Equation (RDE) for a random variable related to N. Unfortunately, our method 
does not give much insights on the heavy-tail phenomena involved in the birth-and- 
assassination process. 

Rumor scotching process 

We now define the rumor scotching process on a graph. It is a nonstandard SIR 
dynamics (see for example [10] or [I] for some background). This process represents the 
dynamics of a rumor /epidemic spreading on the vertices of a graph along its edges. A 
vertex may be unaware of the rumor/susceptible (S), aware of the rumor and spreading 
it as true/infected (I), or aware of the rumor and trying to scotch it/recovered (R). 

More formally, we fix a connected graph G = (V, E) , and let Vy denote the set of 
subsets of V and X = (Vy X {S, I, R}) v . The spread of the rumor is described by a 
Markov process on X. For X = (X v ) ve y G X, with X v = (A v ,s v ), A v is interpreted 
as the set of neighbors of v which can change the opinion of v on the veracity of the 
rumor. If (uv) G E, we define the operations E uv and E v on X by (X + E uv ) w = 
(X - E v ) w = X w ,ifw^v and (X + E uv ) v = (A v U {u}, I), (X - E v ) v = (0, R). Let 
A > be a fixed intensity, the rumor scotching process is the Markov process with 
generator: 

K(X,X + E UV ) = Xl(s u = I)l((u,v)eE)l(s v ^R), 
K(X,X-E V ) = l(s v = I)J2Hsu = R), 

u£A v 

and all other transitions have rate 0. Typically, at time 0, there is non-empty finite set 
of /-vertices and there is a vertex v such that A v contains a i?-vertex. The absorbing 
states of this process are the states without /-vertices. The case when at time 0, A v is 
the set to all neighbors of v is interesting in its own (there, A v does not evolve before 
s v = R). 

If G is the infinite fc-ary tree this process has been analyzed by Kordzakhia [7] 
and it was defined there as the chase-escape model. It is thought as a predator-prey 
dynamics: each vertex may be unoccupied (S), occupied by a prey (I) or occupied by a 
predator (R). The preys spread on unoccupied vertices and predators spread on vertices 
occupied by preys. If G is the Z d -lattice and if there is no /^-vertices, the process is 
the original Richardson's model [13] . With R- vertices, this process is a variant of the 
two-species Richardson model with prey and predators, see for example Haggstrom and 
Pemantle [Bj, Kordzakhia and Lalley |8j. Nothing is apparently known on this process. 

In Section [3], we show that the birth-and-assassination process is the scaling limit, 
as n goes to infinity, of the rumor scotching process when G is the complete graph over 
n vertices and the intensity is X/n (Theorem H]). 
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2 Integral equations for the birth-and-assassination 



process 

2.1 Proof of Theorem [3] for the first moment 

In this paragraph, we prove Theorem [3]^ii). Let X(t) G [0, +oo] be the total number 
of born particles in the process B given that the root cannot die before time t, and 
Y(t) be the total number of born particles given that the root dies at time t. By 

definition, if D is an exponential variable with mean 1 independent of Y, then iV = 

X(0) = Y(D), where the symbol = stands for distributional equality. We notice also 

that the memoryless property of the exponential variable implies X(t) — Y(t + D). 
The recursive structure of the birth-and-assassination process leads to the following 
equality in distribution 

Y(t) £i+j2x(t-&±i+ MM, 

where $ = is a Poisson point process of intensity A and (Xi),i G N, are 

independent copies of X. Note that since all variables are non-negative, there is no 
issue with the case Y(t) = +oo. We obtain the following RDE for the random function 
Y: 

Y(t) = l+ ^Yifa + Di), (3) 

i:£i<t 

where Y i: and Di are independent copies of Y and D respectively. This last RDE is 
the cornerstone of this work. 

Assuming that E\N < oo we first prove that necessarily A G (0, 1/4). For conve- 
nience, we often drop the parameter A in E\ and other objects depending on A. From 
Fubini's theorem, EX(0) = EN = J °° EY(t)e~*cft and therefore EY(t) < oo for almost 
all t > 0. Note however that since 1 1— > Y(t) is monotone for the stochastic domination, 
it implies that EY(t) < oo for all t > 0. The same argument gives the next lemma. 

Lemma 1 Let t > and u > 0, ifE[N u ] < oo then E[Y(t) u ) < oo. 

Now, taking expectation in Q, we get 

EY(t) = 1 + A / / EY(x + s)e~ 8 dsdx. 
Jo Jo 

Let fi(t) = EY(t), it satisfies the integral equation, for all t > 0, 

i*t poo 

f l (t) = l + X e x f^e-'dsdx. (4) 

Jo Jx 

Taking the derivative once and multiplying by e - ', we get: /{(t)e _< = A f°° fi(s)e~ s ds. 
Then, taking the derivative a second time and multiplying by e*: f"{t) — f[(t) = 
—Xfi(t). So, finally, /i solves a linear ordinary differential equation of the second order 

x" - x' + \x = 0, (5) 
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with initial condition x(0) = 1. If A > 1/4 the solutions of (JSJ) are 

x(t) = e t/2 (cos(W4:X - 1) + asin(t v / 4A- 1)), 

for some constant a. Since is necessarily positive, this leads to a contradiction 
and EiV = oo. Assume now that < A < 1/4 and let 

A = Vl -4A , a = and (3 = (6) 

(a, f3) are the roots of the polynomial X 2 — X + A = 0. The solutions of (j3j) are 

x a (t) = (l- a)e at + oef* 
for some constant a. Whereas, for A = 1/4, a = 1/2 and the solutions of (jSJ) are 

z a (t) = (at + l)e' /2 . 

For < A < 1/4, we check easily that the functions x a with a > are the nonnegative 
solutions of the integral equation (J4]). 

It remains to prove that ifO<A<l/4 then EiV < oo and f\(t) = e at . Indeed, 
then EiV = J °° fi(t)e~ t dt = (1 — a)" 1 as stated in Theorem [3](ii). To this end, define 

f[ n \t) = Emm(Y(t),n), from ©, 

min(y(t), n) < st 1 + ^ min^fo + A), n). 

i-4i<t 

Taking expectation, we obtain, for all t > 0, 

ft /»oo 

f[ n \t) < 1 + A / e x f[ n \s)e- s dsdx. (7) 

J0 Jx 

We now state a lemma which will be used multiple times in this paper. We define 

7(A) = (l + A)/(l-A)=/3/a. (8) 

Let 1 < u < 7 (or equivalently A < u(u + l) -2 ), we define TC U , the set of measurable 
functions h : [0, oo) — > [0, oo) such that h is non-decreasing and sup <>0 h(t)e~ uat < oo. 
Let C > 0, we define the mapping from H u to TC U , 

pt poo 

Ce uat + A / e x h(s)e- s dsdx. 

JO Jx 

In order to check that $ is indeed a mapping from 7i u to 7i u , we use the fact that if 
1 < u < 7, then ua < 1. Note also that if 1 < w < 7, then «a — A — u 2 a 2 > 0. If 
A = 1/4, we also define the mapping from TCi to Hi, 

/*£ /*oo 

$:/ii^l + - / e s / h(s)e~ s dsdx. 
4 Jo A 

(recall that for A = 1/4, a = 1/2). 
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Lemma 2 (i) Let 1 < u < 7 and f G 7Y U suc/i i/ioi / < ^(/). T/ien /or a// £ > 0, 

/w < c ^\~^^ - C 

«a — A — u^cr — A — irar 

(nj // A = 1/4 and f E Hi is such that f < $(/), tnen /or a// 1 > 0, 

/(£) < e'/ 2 . 

Before proving Lemma [2j we conclude the proof of Theorem ^ii). For < A < 1/4, 
from ([7]), we may apply Lemma [2]^i) applied to 1 < u < (3/a, C = 1. We get that 

f[ n \t)<C u e aut 

for some C u > 0. The monotone convergence theorem implies that /i(f) = lim^oo f[ (f) 
exists and is bounded by C u e aut . Therefore f\ solves the integral equation (j3J) and is 
equal to x a for some a > 0. From what precedes, we get x a (t) < C u e aut , however, since 
au < (3, the only possibility is a = and fi(t) = e at . 

Similarly, if A = 1/4, from Lemma G9^)> fi(t) < e'/ 2 . This proves that fi is finite, 
and we thus have fi = x a for some a > 0. Again, the only possibility is a = since 
^a(f) < e'/ 2 implies a = 0. 

Proof of Lemma\^ (i). The fixed points of the mapping \1/ are the functions h a ,b such 
that 

ua(l — ua) 
ua — A — -u 2 a 2 

with a + b + C^~^z^5~2 — C- The only fixed point in H u is h* := /i a »,o with a* = 
— CA/ (ua — A — w 2 a 2 ). Let C n denote the set of continuous functions in H u , note that 
$ is also a mapping from C u to C u . Now let go E C„ and for k > 1, = \&(ofe-i). We 
first prove that for all f > , lim^ <7fc(f) = /i*(f). If 1 < u < 7 then wa(l — ua) > A and 
uq-a-"v ^ s P os iti ve - We deduce easily that if g (t) < Le uat then g\{t) = ty(g)(t) < 
c uat + l\ ( uat _ 1) < L x e uat . with Li = (C + n A r ). By recursion, we obtain 

ua(l— ua) \ ) — L ' J- V 1 iitt(l-ua) > J ' 

that limsup fc Ofc(^) < L^e^, with L ro = Cuct(l — uot)/(ua — A — w 2 a 2 ) < oo. From 
Arzela-Ascoli's theorem, (gk)keN is relatively compact in C u and any accumulation point 
converges to h* (since h* is the only fixed point of ^ in C u ). 

Now since / G H u , there exists a constant L > such that for all t > 0, f(t) < 
go(t) := Le uat . The monotonicity of the mapping \1/ implies that \l/(/) < *&{go) = gi- 
By assumption, / < *$?(f) thus by recursion / < lim n o n = h*. 

{%%). The function xo(t) = e t//2 is the only fixed point of $ in Hi. Moreover, if 
g(t) < Ce l l 2 then we also have $(g)(t) < Ce*/ 2 . Then, if g is continuous, arguing as 
above, from Arzela-Ascoli's theorem, ($ fc (g))fc g N converges to Xq. We conclude as in 
(i). " ' □ 

2.2 Proof of Theorem E](i) 

We define f p (t) = E\[Y(t) p ]. As above, we often drop the parameter A in Ea and other 
objects depending on A. 



,(t)=ae^ + be^ + C ^\- Ua 2 ) 2 e^, 
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Lemma 3 Letp > 2, there exists a polynomial Q p with degree p such that for allt > 0, 

(i) IfX G (0,p(p+l)- 2 ), then f p (t) = Q p (e at ). 

(ii) If A > p(p + l)~ 2 , then f p (t) = oo, 

Note that if such polynomial Q p exists then Q p (x) > 1 for all x > 1. Note also that 
A G (0,p(p + 1)~ 2 ) implies that p < 7 = (3/ a (where 7 was defined by ([8])), and thus 
pa < (3 < 1. Hence Lemma [3] implies Theorem [3]^ i) since E[iV p ] = J f p (t)e^ t dt. 

Let k p (X) denote the p th cumulant of a random variable X whose moment gen- 
erating function is defined in a neighborhood of 0: lnEe 6 ^ = ^2 p>0 n p (X)8 p /p\. In 
particular no(X) = 0, ki(X) = EX and k 2 (X) = VarX. Using the exponential for- 
mula 



valid for all non-negative function h and iid variables (Zi),i G N, independent of 
$ = {£i}i £ N a Poisson point process of intensity A, we obtain that for all p > 1, 



Due to this last formula, it will be easier to deal with the cumulant g p {t) = K p (Y(t)). 
By recursion, we will prove the next lemma which implies Lemma [31 

Lemma 4 Let p > 2, there exists a polynomial R p with degree p, positive on [1, 00) 
such that, for all t > 0, 

(1) IfXe (0,p(p+l)- 2 ), then f p {t) < 00 and g p {t) = R p (e at ). 

(ii) If A > p{p + l)~ 2 , then f P it) = 00, 

Proof of Lemma^ In §2.1[ we have computed f p for p = 1 and found i?i(x) = x. Let 
p > 2 and assume now that the statement of the Lemma H] holds for q = 1, ■ ■ ■ ,p — 1. 
We assume first that / p (£) < 00, we shall prove that necessarily A G (0,p(p + 1)~ 2 ) 
and g p (t) = R P ie at ). Without loss of generality we assume that < A < 1/4. From 
Fubini's theorem, using the linearity of cumulants in ([3]) and fflUj) . we get 



(note that Fubini's Theorem implies the existence of f p (s) for all s > 0). From 
Jensen inequality f p (t) > (?i(t) p = c pat and the integral e pas e~ s dsdx is finite if 
and only ii pa < 1. We may thus assume that pa < 1. We now recall the identity: 
EX P = J2 n ri/G7r K \i\{X)> where the sum is over all set partitions of {1, • • • ,p}, I G n 
means I is one of the subsets into which the set is partitioned, and |/| is the cardinal 




(9) 




(10) 




(11) 
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of I. This formula implies that EX P = k p (X) + S p _ 1 (/t 1 (X), • • • ,k p -i(X)), where 
E p _i(xi, • • • , x p _i) is a polynomial in p — 1 variables with non-negative coefficients and 
each of its monomial n£=i satisfies nii£ = p. Using the recurence hypothesis, 
we deduce from ( fTTl) that there exists a polynomial R p (x) = Y7k=i r 'k xk °f degree p 



with r p > such that 



g p (t) = X f Q eX f x ( K 9 P (s)e- s + R p (e as )e- S ^j dsdx 

P \ ft poo 

(recall that pa < 1). Now we take the derivative of this last expression, multiply by e _< 
and take the derivative again. We get that g p is a solution of the differential equation: 

p 

x" -x' + Xx = -^Xr k e ka \ (13) 
fc=i 

with initial condition x(0) = 0. Thus necessarily g p (t) = ae at + be /3t +ip(t) , where ip(t) is 
a particular solution of the differential equation (fT3l . Assume first that A 7^ p(j9+ 1) -2 , 
then it is easy to check that (p + 1)A — pa and p(p + 1)~ 2 — A are different from and 
have the same sign. Looking for a function if of the form ip(t) = EtUi Ck^ hat gives 
c k = -Xr k (k 2 a 2 -ka + A)" 1 = Xr k (k - 1) -1 ((A; + 1)A - ka)' 1 . If A > p(p + l)" 2 then 
pa > (3 and the leading term in g p is c p e pat . However, if A > p{p + 1)~ 2 , c p < and 
thus <7 P (£) < for t large enough. This is a contradiction with Equation (jlip which 
asserts that is positive. 

We now check that if < A < p(p + 1)~ 2 then f p (t) is finite. We define fi n \t) = 
E[min(F(t), n) p }. We use the following identity, 

p— fc— 1 




Then from (j3J) we get, 

(Y(t) - If £ (14) 



E Y ^ + D ^ + E E ( p fc x ) + ( E + ^ 

&<t &<t k=o \ ' \ijHi<t 

The recursion hypothesis implies that there exists a constant C such that fk{t) 
Qk(e at ) < Ce kat for all 1 < k < p - 1. Thus, the identity F(t) p = (Y(t) - l) p 
EG K)(-ir*n*)* g^es 



p— fc— 1 



fcat 



/W(f) < E[min(F(t)-l,^] + ^^Ce 

fc=0 ^ ' 

< E[min(F(t) - 1, n) p ] + de pat . 
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From the recursion hypothesis, if 1 < k < p — 1, 

rt rt poo 

/ E[Y(x + D) k }dx = / e x / f k {s)e- s dsdx = Q k (e at ) < Ce k 

JO JO Jx 



for some constant C > 0. We take the expectation in tHM and use Slyvniak's theorem 
to obtain 



*t 

ti n \t) < C ie pat + \l e x I f ( p n \s)e- s dsdx 

Jo Jx 



+ A /E(% 1 ) E ^( a; +^) fe+1 ] E 

< C ie pat + \ / e x / f ( p n \s)e- s dsdx 

Jo Jx 



p-k-1 



dx 



+ A E ( p k ^www - ir*- 1 ] 

/»£ /»oo 

< C 2 e pai + A / e 1 / f i p n \s)e- s dsdx 

Jo Jx 

So finally for a suitable choice of C, 

ft fOQ 

\f( n \t)<Ce pat + \ e x f { p n \s)e~ s dsdx. 

Jo Jx 



(15) 



From Lemma [21 fp n \t) < C'e pat , and, by the monotone convergence theorem, <7 P (£) < 
/ P (*) < C"e pai . From what precedes: g p (t) = ae at +be pt '+tp(t), with <p(t) = Y7 k= i c k e kat , 
with c p > 0. If 6 > 0, since A > p(p + 1)~ 2 then pa < (3 and the leading term in g p 
is fre' 3 * which is in contradiction with g p (t) < C'e pat . If b < 0, this is a contraction 
with Equation pip which asserts that g p (t) is positive. Therefore 6 = and g p (t) = 
ae at + <p(t) = R p (e at ). 

It remains to check that if A = p(p + 1)~ 2 then for all t > 0, f p (t) = oo. We have 
proved that, for all A < p(p + 1)~ 2 , p (t) = Up(A)(p - l) _1 ((p + 1)A - pa)~ l e pat + 
S p -i(e at ), where S p -i is a polynomial of degree at most p—1 and u p (X) > 0. Note that 
limx-|- p (p+i)-2 (p + 1)A — pec = 0. A closer look at the recursion shows also that u p (X) is 
a sum of products of terms in A and X(£ — 1) _1 ((£ + 1)A — ia)^ 1 , with 2 < £ < p — 1. 
In particular, we deduce that lim A -|- p ( p+1 )-2 w p (A) > 0. Similarly, the coefficients of S p _i 
are equal to sums of products of integers and terms in A and X(£ — 1)~ 1 ((£+1)A — ia) -1 , 
with 2 < £ < p — 1. Thus they stay bounded as A goes to p(p + 1)~ 2 and we obtain, 
for all t > 0, 

liminf f„(t) > lim o„(t) = oo. (16) 

Now, for all t > 0, the random variable Y(t) is stochastically non- decreasing with A. 
Therefore EA[V(i) p ] is non-decreasing and (|T6|) implies that E 1 / 4 [y(t) p ] = oo. The 
proof of the recursion is complete. 

□ 



10 



2.3 Proof of Theorem E](iii) 

In this paragraph, we prove Theorem EJ^iii). Let A G (0, 2/9), recall that f2{t) = EY(t) 2 
and g2(t) = Var(F(i)). From (fTTj) applied to p — 2, 

02 (t) = A / / # 2 (£ + s)e~ s + / 2 (x + s)e~ s dsdx. 
Jo Jo 

Since = e at and a 2 = a — A, #2 satisfies the integral equation: 



\ /*£ /*00 

We deduce that g 2 solves an ordinary differential equation: 

x" - x' + Xx = -Xe 2at , 
with initial condition x(0) = 0. Thus #2 is of the form: g2(t) = ae at + be^ + 



A n 2at 



3X-2a 

with a + b + 3X ^ 2a = 0. From Lemma HI b = so finally 



g J t ) = (e 2at -e at ) and f 2 (t) = 2 2A - e 2 ^ e at . 

y w 3A-2a V ; J W 3A-2a 3A - 2a 

We conclude by computing EA^ 2 = J e~ t f 2 (t)dt. 
2.4 Proof of Theorem H 



As usual we drop the parameter A in E^. From (jSJ), we have 7(A) = 1 2A ^ 1 — ■ To 
prove Theorem [21 we shall prove two statements 

If E[N U ] < 00 then u < 7, (17) 
If 1 < u < 7 then ELY"] < 00. (18) 

2.4.1 Proof of (TIT]). 

Let m > 1, we assume that E[A^ U ] < 00. From Lemma [Hand ([3]), we get 



E[y (*n = e (1 + j^y^+Di)) 



Let f u {t) = E[Y(t) u ]. Taking expectation and using the inequality (x + y) u > x u + y u , 
for all positive x and y, we get: 



/«(*) > 1 + A / E/„(x + D)cfc 
Jo 



n t POO 

> 1 + A / e x f u (s)e- s dsdx. (19) 

</x 
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From Jensen's Inequality, f u (t) > }\{t) u = e uat . Note that the integral e aus e~ s ds 
is finite if and only if u < a" 1 . Suppose now that 7 < u < a -1 . We use the fact: if 
u > 7 then u 2 a 2 — ua + A > 0, to deduce that there exists < e < A such that 

u 2 a 2 -ua + X> e. (20) 

Let A = A — e, a = a (A), /3 = /3(A), we may assume that e is small enough to ensure 
also that 

ua > (3. (21) 

(Indeed, for all A G (0, 1/4), «(A)t"(A) = /3(A) and the mapping A 1— > /3(A) is obviously 
continuous). We compute a lower bound from (fill as follows: 

/*£ /*oo /»£ ^00 

> 1 + A / e x / f u (s)e~ s dsdx + e e x f u (s)e~ s dsdx 

Jo Jx Jo Jx 

t*t /*oo /»£ /»oo 

> 1 + A / e* / f u (s)e- s dsdx + e / e x / e uas e~ s dsdx 

Jo Jx Jo Jx 



pt poo 

> l + C(e uat + \ e x f u (s)e~ s dsdx, (22) 

Jo Jx 

with C = e(ua(l — «a)) _1 > 0. We consider the mapping : h 1— > 1 + C , (e Mat — 
1) + A J* e x h(s)e~ s dsdx. \1/ is monotone: if for all t > 0, /ii(t) > h 2 {t) then 
for all t > 0,^(/ii)(t) > V(h 2 )(t). Since, for all t > 0, f u (t) > *(/„)(*) > 1, we 
deduce by iteration that there exists a function h such that = ty(h) > 1. Solving 
/i = ^(h) is simple, taking twice the derivative, we get, h" — h' + \h = —ee pat . Therefore, 
h = ae at + be 131 — e(u 2 a 2 — ua + \)~ l e uat for some constant a and b. From ( !2~TI) the 
leading term as t goes to infinity is equal to — e(u 2 a 2 — ua + \)~ 1 e uat . However from 
( |20|) . — e(w 2 a 2 — ua + A) -1 < and it contradicts the assumption that h(t) > 1 for all 
t > 0. Therefore we have proved that u < 7. 

2.4.2 Proof of (TI8]). 

Let f [ u\t) = E[ min(F(t), n) u ], we have the following lemma. 
Lemma 5 There exists a constant C > snc/i that for all t > 0: 

fi n) (t) < Ce uat + A / e* / fi n \s)e- s dsdx. 



JO Jx 

The statement ffl8l) is a direct consequence of Lemmas [2] and [5j Indeed, note that 
fi n) < n u , thus by Lemma El for all t > 0, f [ u\t) < de uat for some positive constant 
C\ independent of n. From the Monotone Convergence Theorem, we deduce that, for 
all t > 0, f u {t) < C\e uat . It remains to prove Lemma [51 

Proof of Lemma The lemma is already proved if u is an integer in (Tl5l) . The general 
case is a slight extension of the same argument. We write u = p — 1 + v with v G (0, 1) 
and p G N*. We use the inequality, for all yi > 0, 1 < i < N, 

. i=i / j=i fc=o v 7 Vjyi 
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(which follows from the inequality (J^ yj) v < ^2 y%). Then from ([3]) we get the stochas- 
tic domination 

(Y(t)-l) u < st J^Yifa + DiT 

From Lemma [3j there exists C such that for all 1 < k < p — 1, /&(£) < Ce kat and and 
J E[F(x + D) fe ](ix < Ce kat . Note also, by Jensen inequality, that for all 1 < A; < p — 2, 
fk+v{t) < / p -i(t) (fc+,,)/(p_1) < (7e (fc+,,)Q! *. The same argument (with p replaced by u) 
which led to ( 1T51) in the proof of Lemma H] leads to the result. □ 



2.5 Some comments on the birth-and-assassination process 
2.5.1 Computation of higher moments 

It is probably hard to derive an expression for all moments of N, even if in the proof of 
Lemma HI we have built an expression of the cumulants of Y(t) by recursion. However, 
exact formulas become quickly very complicated. The third moment, computed by 
hand, gives 

3A - a 3at A(2A - a) 2at ( A(2A - a) 3A - a \ t 

Since N = Y(D), we obtain, 

EN 3 = 6 ( 3A -^) a 6 \(2\-a)a | 1 



(4A-3a)(l -a-3A) (3A - 2a) 2 (l - a - 2A) I- a 
2.5.2 Integral equation of the Laplace transform 

It is also possible to derive an integral equation for the Laplace transform of Y(t): 
Lg(t) = Eexp(— 9Y(t)), with 6 > 0. Indeed, using RDE (j3J) and the exponential 
formula ©, 

L 6 (t) = e~ e exp ^A (EL e (x + D) — l)dx 



e e exp^A J e x J (Lg(s) - l)e s rfs<ix^ 



Taking twice the derivative, we deduce that, for all 6 > 0, Lg solves the differential 
equation: 

x"x — x' 2 — x'x + \x 2 {x — 1) = 0. 
We have not been able to use fruitfully this non-linear differential equation. 
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2.5.3 Probability of extinction 



If A > 1/4 from Corollary [TJ the probability of extinction of B is strictly less than 1. It 
would be very interesting to have an asymptotic formula for this probability as A get 
close to 1/4 and compare it with the Galton- Watson process. To this end, we define 
Tr(t) as the probability of extinction of B given than the root cannot die before t. With 
the notation of Equation Q, 7r(i) satisfies 

tt(*) = E J] n{t + D-ti)=E J] tt(&), 

i:S,i<t+D i:^<t+D 

Using the exponential formula (jUj), we find that the function 7r solves the integral 
equation: 

7r(t) = e* / exp I —(A + l)s + A / n{x)dx I ds. 



After a quick calculation, we deduce that tt is solution of the second order non-linear 
differential equation 

Unfortunately, we have not been able to get any result on the function ir(t) from this 
differential equation. 



3 Rumor scotching in a complete network 
3.1 Definition and result 




We consider the rumor scotching process on the graph G n on {0, ■ • • , n} obtained 
by adding on the complete graph on {1, • • • , n} the edge (0, 1), see FigureEJ Let V n be 
the set of subsets of {0, • • • , n}. With the notation in introduction, the rumor scotching 
process on G n is the Markov process on X n = (V n x {S, I, R}) n with generator, for 
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X — (A4, Sj)o<j< n , 



K(X, X + Eij) = Xn- 1 !^ = I)l( Sj ^ R), 
K(X,X- Ej) = l( Sj = I) \J2 1(* e A,-) j , 

and all other transitions have rate 0. At time 0, the initial state is X(0) = (Xj(0))o<j<„ 
with X (0) = (0, R), Xx(0) = ({0}, I) and for i > 2, X { (0) = (0, S). 

With this initial condition, the process describes the propagation of a rumor started 
from vertex 1 at time 0. After an exponential time, vertex 1 learns that the rumor is 
false and starts to scotch the rumor to the vertices it had previously informed. This 
process is a Markov process on a finite set with as absorbing states, all states without 
/-vertices. We define N n as the total number of recovered vertices when the process 
stops evolving. We also define Y n (t) as the distribution N n given that vertex 1 is 
recovered at time t. We have the following 

Theorem 4 (i) If < A < 1/4 and t > 0, as n goes to infinity, N n and Y n (t) 
converge weakly respectively to N and Y(t) in the birth-and-assassination process 
of intensity A . 

(ii) If X > 1/4, there exists 5 > such that 

liminf P A (iV n > 8n) > 0. 

n 

The proof of Theorem H] relies on the convergence of the rumor scotching process 
to the birth-and-assassination process, exactly as the classical SIR dynamics converges 
to a branching process as the size of the population goes to infinity. 



3.2 Proof of Theorem H 
3.2.1 Proof of Theorem 

The proof of Theorem H] relies on an explicit contruction of the rumor scotching process. 
Let < i < j < n, be a collection of independent exponential variables with 

parameter Xn^ 1 and, for all 1 < i < j, let be an independent exponential variable 
with parameter 1. We set = and = Qj . A network being a graph 
with marks attached on edges, we define K. n as the network on the complete graph of 
{1, • • • ,n} where the mark attached on the edge (ij) is the pair (^j\Dij). Now, the 

rumor scotching process is built on the network /C„ by setting as the time for the 
infected particle % to infect the particle j and Dij as the time for the recovered particle 
% to recover the particle j that it had previously infected. 

The network K, n has a local weak limit as n goes to infinity (see Aldous and Steele 
[3J for a definition of the local weak convergence). This limit network of JC n is /C, the 
Poisson weighted infinite tree (PWIT) which is described as follows. The root vertex, 
say 0, has an infinite number of children indexed by integers. The marks associated to 
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the edges from the root to the children are -Dj)i>i where {£i}j>i is the realization 
of a Poisson process of intensity A on IR + and (A)i>i is a sequence of independent 
exponential variables with parameter 1. Now recursively, for each vertex % > 1 we 
associate an infinite number of children denoted by (i, 1), (i, 2), • ■ • and the marks on the 
edges from i to its children are obtained from the realization of an independent Poisson 
process of intensity A on R + and a sequence of independent exponential variables with 
parameter 1. This procedure is continued for all generations. Theorem 4.1 in |3J implies 
the local weak convergence of JC n to /C (for a proof see Section 3 in Aldous [T]). 

Now notice that the birth-and-assissination process is the rumor scotching process 
on /C with initial condition: all vertices susceptible apart from the root which is infected 
and will be restored after an exponential time with mean 1. 

For s > and £ e N, let JC n [s,£) be the network spanned by the set of vertices 
j G {1, • • • , n} such that there exists a sequence (ii, • • • , ik) with i 1 — 1, if. = j, k < £ 
and max((| i n i ' 2 , • • • < s. If r n is the time elapsed before an absorbing state is 

reached, we get that 1(t„ < s)l(N n < £) is measurable with respect to fC n [s, £]. From 
Theorem 4.1 in [3], we deduce that l(r n < s)l(N n < £) converges in distribution to 
l(r < s)l(N < t) where r is the time elapsed before all particles die in the birth- 
and-assassination process. If < A < 1/4, r is almost surely finite and we deduce the 
statement (i). 

3.2.2 Proof of Theorem [Kii) 

In order to prove part (ii) we couple the birth-and-assassination process and the rumor 
scotching process. We use the above notation and build the rumor scotching process 
on the network JC n . If X = ((A, Si)o<i<n) G X n , we define I(X) = {1 < i < n : s, = /} 
and S(X) = {l<i<n: Si = S}. 

Let X = X n {u) G X n be the state of the rumor scotching process at time u > 0. Let 
i G I(X), we reorder the variables (£,ij^)j£S{x) in non-decreasing order: < • • • < 
Ci^Ljni" -Define £^ = 0, from the memoryless property of the exponential variable, for 
1 < k < \S(X)\, — is an exponential variable with parameter A(|S'(X)| — k + 
l)/n independent of (Q^ — ^ , £ < k). Therefore, for all 1 < k < \S{X)\, the vector 

■ • • is stochastically dominated component-wise by the vector (^i, • • • , 

where {CiL>i is a Poisson process of intensity A(|S'(X)| — k + l)/n on R + (i.e. for all 
< ti < • • • < t k , P(^ } > *i, • • • , > 4) < P(^i >*!,•••,&> 4))- In particular 
if |S(X)| > (l-5)n, with < 5 < 1/2, then (^ n) , • • • , is stochastically dominated 

component- wise by the first \n5\ arrival times of a Poisson process of intensity A(l— 28). 

Now, let 5 > such that A' = A(l - 25) > 1/4. We define S { u\l { u\^u \ as the 
number of S, I, i?-particles at time u > in JC n , and I' u as the number of particles 
"at risk" at time u in the birth-and-assassination process with intensity A'. Let r n = 
inf{ii > : < (1 — 5)n}. Note that if < u < r n then any /-particle has infected 
less than [8n\ S-particles. From what precedes, we get 

S^l(u<r n ) < st n-l' u . 
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So that < st max(n — V u , (1 — 8)n). In particular, since N n > sup u>0 (n — S^), we 
get 

P A (iV n > 8n) > P A /(IimsupJ£ = oo). 

u— >oo 

Finally, it is proved in |2] that if A' > 1/4 then PA^limsup^oo I' u = oo) > 0. 
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